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We carry out a comprehensive theoretical and experimental study of charge injection in Poly(3- 
hexylthiophene) (P3HT) to determine the most likely scenario for metal-insulator transition in this 
system. We calculate the optical absorption frequencies corresponding to a polaron and a bipolaron 
lattice in P3HT. We also analyze the electronic excitations for three possible scenarios under which 
a first- or a second-order metal-insulator transition can occur in doped P3HT. These theoretical 
scenarios are compared with data from infrared absorption spectroscopy on P3HT thin film field- 
effect transistors (FET). Our measurements and theoretical predictions suggest that charge- induced 
localized states in P3HT FETs are bipolarons and that the highest doping level achieved in our 
experiments approaches that required for a first-order metal-insulator transition. 



I. INTRODUCTION 

Upon doping the conductivity of conjugated organic 
polymers can be controlled over an extraordinarily large 
range from insulator to metali Therefore, these mate- 
rials hold great promise for high-performance organic 
electronics^,^ The organic field-effect transistor is a 
benchmark system for exploring the fundamental prop- 
erties of the organic semiconductors as well as for ex- 
ploiting their novel functionalities^ Recently, particular 
interest has been devoted to FETs based on highly regio- 
regular (RR) P3HT (see Fig. |T) &L2&mL dl dZlldZ. d£ 
The charge carrier mobility in these devices has been 
found as high as 0.7 cm 2 /Vs, enhanced by orders of mag- 
nitude compared to that of the more disordered P3HT 
films. There has also been a strong interest in using the 
voltage in the organic FETs to induce a sufficiently high 
density of carriers so that metallic states can be produced 
in the polymers i 15 i 16 

Unlike traditional semiconductors, charge injection 
into organic polymers induces structural deformations in 
the polymer chains and gives rise to self-localized excita- 
tion states inside the band gap. The Peierls instability 
in trans-polyacetylene (PA) leads to solitonsjil or kink- 
like excitations with a two-fold-degenerate ground state. 
However, most of the conjugated polymers are analo- 
gous to cis-PA, which consists of two non-degenerate 
ground states (NDGS); i.e., the two conformations re- 
sulting from the mutation of a single and a double C-C 
bond are not energetically equivalent. Two possible ex- 
citations in NDGS conjugated polymers are polarons or 
bipolarons that are associated to a pair of bound states in 
the band gap split off from the continuum. 17 Experimen- 
tal evidence of similar signatures have been reported in 
P3HT and have been interpreted in terms of both types 
of excitations. Earlier work on chemical doping and pho- 
toinduced absorption in P3HT have assigned the two lo- 
calized intragap states observed to bipolarons,— whereas 
a more recent study of RR-Poly(3-alkylthiophene) con- 
cluded that while chemical doping in P3HT films results 



in localized polarons, photoexcitations give rise to de- 
localized polarons. 9 A charge modulation spectroscopy 
experiment in regio-random P3HT in a metal-insulator- 
semiconductor diode has identified five sub-gap transi- 
tion features, out of which two are associated to opti- 
cal transitions of doubly-charged bipolarons and three 
to singly-charged polarons, suggesting coexistence of the 
two charge species*^ More recently, however, the same 
technique applied to RR P3HT identified only three sub- 
gap features, which were assigned to polarons £211 

Metallic states attested by a finite density of states 
at the Fermi energy have been reported in chemically- 
doped P3HT polymers using photoemission spectroscopy 
and have been attributed to the polaron latticed On 
the other hand, there exists a model which supports the 
bipolaron lattice in NDGS polymers, arguing that the 
metallic states are a consequence of the merging of the 
bipolaron band with either the valence or the conduction 
band. 22 However, in neither of these studies a prediction 
of the critical doping concentration under which P3HT 
becomes metallic has been given. 

Here we report a joint theoretical and experimental 
study of charge-induced electronic excitations in P3HT. 
We present models for three scenarios under which 
charge-induced transitions between the bipolaron lattice 
and the polaron lattice or the merging between the local- 
ized states and the extended states can take place. We 
predict the critical doping thresholds above which metal- 
lic states can be expected in charge-injected P3HT. To 
this end, we apply the continuum Hamiltonian model in- 
troduced by Brazovskii an Kirova 23 (BK) for that takes 
into account the electronic and lattice couplings. We car- 
ried out infrared spectroscopic studies in P3HT- based 
FETs. From these measurements we extract the charge- 
injected doping levels and compare with the predicted 
critical values, which provide a lower bound for metallic 
behavior. The analysis suggests that the charge carriers 
in our P3HT FET are bipolarons and the highest charge 
injection level reaches the proximity of the critical dop- 
ing value that is needed for metallic states in P3HT. We 
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FIG. 1: Schematic diagram of RR P3HT. 

also discuss our results in view of very recent experimen- 
tal evidence of field-induced metallic behavior in P3HT 
FETs,l— and a structurally similar polymer with poly- 
electrolyte gate dielectric— 

The paper is organized as follows. In Sec. [ill we briefly 
review the solution of the BK model as derived in the lit- 
erature and provide an analytical formulation for band- 
widths of polarons and bipolarons. In Sec. IIII| the solu- 
tions are applied to three scenarios involving polaron and 
bipolaron states in P3HT. From this study the critical 
dopings required to achieve metallic states in P3HT are 
predicted. In Sec. HV| we present experimental results on 
charge injection in P3HT-based transistors and we com- 
pare them with the theoretical predictions in Sec. [V] 

II. THEORETICAL FORMULATION 

A. Polarons and bipolarons 

The Su-Schrieffer-Heeger (SSH) 25 electron-phonon 
coupling model and its continuum version 2 ^ have been 
widely used to study structural excitations (solitons, po- 
larons) in degenerate polymers such as trans-PA. In 
NDGS polymers, however, the degeneracy is broken. 
A generalization of the SSH model by Brazovskii and 
Kirova has been proposed to take into account the one- 
electron energy gap in the absence of dimerization. The 
solution of this model has been derived and discussed 
elsewhere ) 23 ^ 27 ' 28 here we briefly summarize the main re- 
sults. 

The gap parameter contains both an intrinsic and an 
extrinsic contribution A(x) = A e + Ai(x), where Ai(x) is 
associated with the dimerization due to the 7r-electron- 
lattice interaction and A e is a constant gap parameter 
that arises from the er-bond orbitals. The BK continuum 
Hamiltonian is written as 



H = 



dx 



ttXvf 
dx^\x) 



u . 

~ WF<J ' i 'dx + ( x > ai 



(1) 



Here, the first term describes the interaction of the elec- 
trons with the lattice deformation Aj(x), ^>(x) is the two- 
components fermion field, A is the dimensionless electron- 
phonon coupling strength, and Dp is the Fermi velocity 
of the 7r-electrons. The second term describes the self- 
energy of the 7r-electrons in a field of a dimerization po- 
tential A (a;), while <Xj are the Pauli matrices. 

In the continuum limit, the analytical solution of the 
BK Hamiltonian represents a polaron or a bipolaron, 

A(x) = A + kvF{ta,nh{k(x — x )]—ta,nh{k(x+xo)]}, (2) 

with 

kvF 



tanh(2fciEo) 
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where W is the width of the 7r-band and 7 = is the 
confinement parameter that describes the extent of the 
energy degeneracy breaking. 

There are two bound states corresponding to the po- 
laron solution in Eq. ([2|) with eigenvalues symmetrically 
placed in the gap at ±wo where 



u> 



(5) 



and £0 = a 5 - is the parameter that describes the spa- 
tial extent of the soliton kink. Minimizing the energy 
of Eq. (TT]) with respect to wo yields the following self- 
consistent condition for k: 



kv F . -1,^0, . 7T 

— 7 = sin (tH + t^h 



(6) 



where n± are the occupation numbers of the discrete 
energy levels. For a singly-charged system such as an 
electron polaron or a hole polaron, n_|_ — n_ = — 1, 
and for a doubly charged system such as a bipolaron, 
n + — n_ = 0. The polaron (bipolaron) formation energy 
can be obtained from the total energy given in the above 
configuration 28 
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B. Bandwidth of polarons and bipolarons 



As the doping level increases, the charges stored on 
individual polaron states overlap and give rise to the po- 
laron band. The bandwidth can be calculated from the 
off-diagonal element of the BK Hamiltonian between un- 
perturbed states on each polaron^ Here we follow the 
approach of Onodera3£ in which the interaction between 
two hole-doped polarons approaching each other are con- 
sidered. When the doping level increases, the overlap of 
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electron wave functions trapped in the two polaron dips 
in A(x) gets stronger and the two bound state energy lev- 
els ±ojQ are each split into two levels ±wi and ±ui2 (bond- 
ing and anti-bonding states). The system can lower its 
energy by accommodating both holes on the lowest level. 

In a tight -binding approximation, the bandwidth of 
the polarons can then be calculated from the difference 
between the two discrete levels. The gap parameter for 
the two polaron or two bipolarons is found to be^ 



A(x) = A + 



k\ coth(fcia; — f3\) — k 2 tanh(/c2^ — P%) 
(k\ - k\)v F 



k\ cath-ikxx + (3) — k 2 tanh(fc2X + (3 2 



(8) 



where 



tanh2/3 J = fc^ (i = l,2), (9) 
and k\ and k 2 are related to the discrete energy levels by 



(10) 



As the two polarons approach each other, k\ — k 2 be- 
comes nonzero and is asymptotically related to the sep- 
aration d between the two polarons via 



fci 



4fce 



-<lk 



(11) 



at large distances. Therefore the separation between the 
bonding and antibonding states can be calculated from 
Eq. (I10|) which now becomes 



(fcl 2 - kl)v 2 F + (col - u\) = 0. 



(12) 



If the separation d is large, ki — k% is small, we can rewrite 
this relation as 



k(k\ — k 2 )v F + u>o(lui — w 2 ) = 0. 



(13) 



Replacing d by the dopant concentration y = a/d, where 
a is the projected bond length along the chain axis, we 
obtain the bandwidth of the polarons 



W p (y) 



. 4 ( Wo - W2 )^8^ e -^, 

UJQ 



(14) 



In a similar fashion, the two-bipolaron system has been 
solved,— and an asymptotic expression of the bipolaron 
bandwidth similar to that for the polarons can be ob- 
tained, 



7T V Wn / Wn 



(15) 



Expressions (|14[) and (|15|) represent our starting point 
for the discussion that follows. 
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FIG. 2: Schematic diagram of the localized gap states in a 
P3HT polymer with p-type doping showing three scenarios 
for the metal-insulator transition including: a) the crossover 
between the bipolaron and the polaron formation energy (the 
lower polaron band is partially filled) upon doping; b) the 
polaron bands broaden to merge with the continuum; c) the 
bipolaron band broaden to merge with the continuum. The 
shades are indicative of occupied electronic bands. 



III. THEORETICAL RESULTS 

In NDGS polymers, doped or injected charges can lead 
to the formation of either polarons or bipolarons, de- 
pending on the formation energy of the excitations and 
the relative energy difference between the two configura- 
tions upon doping. Here we apply the solutions to the BK 
model to predict the effect of the doping concentration on 
the electronic excitations in P3HT. We considered three 
types of band diagrams as shown in Fig. O They corre- 
spond to a) a first-order transition between a bipolaron 
and a polaron lattice, b) a gap closure of the polaron band 
and the conduction or the valence band, and c) a gap clo- 
sure between the bipolaron band and the conduction or 
the valence band. In the latter two cases, we look for the 
condition that the polaron (bipolaron) bands merge with 
the extended states and we expect these transitions to be 
of second order. 



A. Optical absorption energies 

Before proceeding to calculate the critical doping con- 
centration that can give rise to metallic states, we calcu- 
late the localized energy levels and the formation energy 
for the polarons and bipolarons of P3HT polymers. The 
two excitation states correspond to a pair of symmetri- 
cally localized levels^ in the gap that can be calculated 
from Eq. ([S]). To solve for k, we apply the self-consistency 
condition Eq. ([6]), 
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TABLE I: The Eigenvalues, the subgap optical transition en- 
ergies and the formation energies of a polaron and a bipolaron 
in P3HT calculated from the BK model. The values in the 
parenthesis are calculated using 7 = 0.3. The units are eV. 



excitations 




0)1 


0J2 


0)3 


E 


Polaron 
Bipolaron 


0.89 (0.87) 
0.59 (0.52) 


0.16 
0.46 


1.94 
1.74 


1.78 


1.0 
1.64 



ergy of two single polarons, as also verified by the for- 
mation energies that we have calculated in the preceding 
subsection. Therefore, upon doping, the charges will be 
stored in bipolarons rather than the polarons. However, 
due to the partial filling, the energy per polaron changes 
with doping, whereas the bipolaron energies remain un- 
changed because the bipolaron bands are either doublely 
occupied or empty. The condition to be satisfied by this 
transition is 



and 



1 - z 2 



(17) 



for the polaron and the bipolaron configuration, respec- 
tively, where z — fc£o- The confinement parameter 7 
can be either estimated from the energy splitting be- 
tween the doping-induced gap states^ or by solving for 
the phonon propagator from Raman resonance dispersion 
data^ In either method, 7 is estimated to be 0.3 ~ 0.4 
in P3HT (we discuss later the effect of disorder on this 
parameter)^ Solving the above two equations, we find 
that z p is in the range 0.553 — 0.517 and Zb in the range 
0.867-0.829. 

The energy gap of P3HT has been determined to be 
2Ao = 2.1 eV.~ Together with the above values of z, 
we find the localized energy levels too that correspond 
to the polaron and the bipolaron. The singly-occupied 
polaron states allow three optical absorption transitions, 
namely, a transition between the valence band and the 
— loq level (u>x), a transition between the valence band 
and the +wo level (u^) and a transition between the two 
localized levels (w 3 )i21 For the bipolarons instead, the 
last type of transition between the two localized levels 
is forbidden because the localized states are either both 
empty or doubly occupied. Using the value of 7 = 0.4 
in P3HT, we calculated the optical transition frequencies 
corresponding to the three transitions for the polaron and 
the two allowed transitions for the bipolarons in Table U 

Substituting the self-consistent condition into Eq. J?}, 
we also find that the formation energy of a polaron and 
a bipolaron can be expressed as E P:( j = —AqN(z) with 
N(z) = z + 7tanh~ 1 (z) and their values are shown in 
Table H 



B. Transition between bipolarons and polarons 

We here discuss the transition between the bipolaron 
lattice and the polaron lattice (see Fig. [2^). This first- 
order phase transition in NDGS conjugated polymers has 
been previously proposed without, however, an explicit 
calculation of the critical doping density 29 i 31 This tran- 
sition is similar to the energy crossover between the soli- 
ton lattice and the metallic polaron lattice, as proposed 
for trans-PA by Kivelson and Heeger^ The formation 
energy of a bipolaron is smaller than the formation en- 



W p (y) = tt(2E p - E b ), 



(18) 



where W p (y) is the bandwidth of the polarons we have 
derived in Section [TTJ Substituting the formation energy 
of a polaron and a bipolaron in the above equation, we 
get 



W p (y) = 4A [2N(z p ) - N{z b )}. 



(19) 



Using the polaron bandwidth, Eq. (jT4|) , we can solve for 
the critical dopant concentration 



l-z%[2N(z p )-N(z b )] 



where £0 is the characteristic length that can be deter- 



mined from £0 



w „ 29 

2A, 



a^sA with a the lattice spacing and 



W the 7T band width. Since the tt bandwidth is a rather 
universal constant for most of conjugated polymers, we 
have taken a typical value W ~ 10 eV— and obtained 
£0 ~ 4.7a. Therefore we find that y c ~ 0.092. 



C. Polaron gap closure 

A second possible scenario corresponds to the case in 
which the polaron energy levels stay unchanged upon 
doping but the bandwidth of the polaron states broadens 
(see Fig. (2b)- When the doping is sufficiently large, the 
polaron band can merge with the extended states. This 
situation can be realized as long as the condition 



W p (y) = 2uj x 



(21) 



is satisfied. Substituting this in the expression for the 
polaron bandwidth we obtain 



l-zj 



-Ww = 2lui 



(22) 



and thus 



£0 



4A z 2 



(23) 



Using the lowest absorption energy of a polaron state 
found in Sec. IIII A[ i.e., oj\ = 0.16 eV, we find that the 
critical doping concentration is y c — 0.051. 
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D. Bipolaron gap closure 



Similar to the case of the polaron, merging between 
the bipolaron band and the valence band (see Fig. [2p) 
can also occur if the condition 



W b (y) = 2uj! 



(24) 



is satisfied. From the bipolaron bandwidth Eq. (| 1 5 [) . the 
above condition becomes 

7T 1-Zt \ I- 4 J 



hence 



2z b a f 7twi(1 - zlf 
Vc = — { hi ■ 



£o I 16A z 3 (l-z 2 + 7 ) 



(26) 



Using the result uj\ = 0.46 eV and the value of z b for a 
bipolaron in Sec. IIII A[ we find that the critical dopant 
concentration for this transition is y c = 0.091. 



IV. INFRARED SPECTROSCOPIC RESULTS 

In order to study the insulator-metal transition in 
conjugated polymers in functional organic FET devices, 
we carried out systematic investigations of P3HT based 
FETs employing Infrared (IR) spectroscopy. The critical 
advantage of IR spectroscopy is that it directly probes the 
electronic excitations associated with the injected carri- 
ers in organic FET devices^ Combined with theoreti- 
cal analysis in Sec. IIII1 these experiments uncover the 
nature of the electronic excitations in P3HT. Moreover, 
the carrier density in the accumulation layer of P3HT 
can be evaluated from the analysis of IR spectroscopic 
data, which provides new insights into the insulator- 
metal transition in polymers. 

Large area organic FET device a 33 i 34 (> 1 cm 2 ) with 
gate insulator deposited on n-Si were investigated in 
this work. In these devices, source and drain Au elec- 
trodes (with a spacing of 50-200 fim) were patterned 
on a 200-nm-thick SiC>2 gate insulator followed by the 
deposition of a 4-6 nm-thick P3HT film. The inset of 
Fig. El displays a schematic of the cross-section of the de- 
vice. These devices has transport mobility of about 0.18 
cm 2 V _1 s _1 and breakdown voltage in excess of — 100V. 35 
In these bottom-contact FET devices, an applied gate 
voltage induces an accumulation layer^& in P3HT as 
well as in n-Si. The former accumulation layer con- 
stitutes the p-type conducting channel between the Au 
electrodes. This channel is not obscured by any other 
interfaces and is therefore well suited for the spectro- 
scopic studies of the accumulation layer in the polymer 
film from far-IR to near-IR with the latter cut-off im- 
posed by the band gap of the Si substrate^ 3 - We stud- 
ied changes of transmission as a function of applied gate 
voltage Vgs normalized by the transmission at Vgs=0: 



100 
<~n — 



co (cm" 1 ) 



1000 
<-n 



7000 



Poly(3-hes^lfh.iophene) (~4 nm) 
Au Au Au Au 
Gateinsulator 

n-doped Si (0.5 mm) 




7000 



FIG. 3: (Color online) Top panel: the 1-T(Vgs)/T(0V) spec- 
tra for a representative Si02 based organic FET device under 
applied gate voltages Vgs- Solid curves: experimental spec- 
tra at several gate voltages. Blue dashed curve: the spectrum 
at — 100V obtained from fitting the low-energy absorption of 
the accumulation layer in n-Si using the DL model as de- 
scribed in the text. Top inset: a schematic of the cross-section 
of the FET devices. Bottom panel: The mid-IR voltage- 
induced absorption spectra Aad for the P3HT layer, which 
are obtained by subtracting absorption due to the accumula- 
tion layer in n-Si from the absorption spectra of the device as 
detailed in the text. 



T(w,Vgs)/T(w,Vgs=0). Note that in the same sample 
structure without a P3HT layer, we did not observe any 
gate voltage induced changes in the IR spectrum. This 
suggests that no accumulation layer forms in the region 
between the gold grid electrodes without the P3HT layer 
on the top of the insulator. The source and drain volt- 
age is zero in most measurements. All the data reported 
here were recorded at room temperature with a spectral 
resolution of 4 cm -1 . 

We start by analyzing the absorption spectra Aad = 
1-T(V GS )/T(0V) of the organic FET devices displayed 
in Fig. [3] Here Aa is the change of the absorption coeffi- 
cient of P3HT and d is the thickness of the accumulation 
layer where absorption is formed. These spectra show 
several voltage-induced features: i) a strong absorption 



in far-IR frequencies (lj < 500 cm 



sharp reso- 



nances in the 1000 - 1500 cm 1 region, and iii) a broad 
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TABLE II: The parameters of the Drude-Lorentz fitting 
(Eq. (|27p ) of the absorption of the accumulation layer in n-Si 
of the FET device at V G s = -100V. We use the "effective 
2D" oscillator strengths u> PD and u>f> associated with 2D car- 
rier densities N/ and Nm, respectively. 



2 



2.57cm" 1 
0.74 x 10 13 . 



'cm 



132cm" 1 
0.96cm" 1 



Nm 

r 

Wo 



0.28 x 10 13 cm" 2 
1400cm" 1 
360cm" 1 
10 



band centered around 3500 cm" 1 . A gradual develop- 
ment of these features with increasing gate voltages Vgs 
suggests that i-iii) are intimately related to the forma- 
tion of charge accumulation layers on both sides of the 
oxide. Features ii) and iii) are spectroscopic signatures 
of electrostatically doped carriers in P3HT. 33 Sharp res- 
onances in the 1000 - 1500 cm _1 range result from the 
IR active vibrational modes (IRAVs)^ i.e Raman modes 
made IR active by distortions of the polymer backbone 
caused by the self-localized charges. The broad absorp- 
tion band centered around 3500 cm" 1 is usually ascribed 
to polaron or bipolaron states positioned in the gapii; 
the lineshape of the broad band does not allow one to 
discriminate between the two types of excitations. Both 
the IRAV modes and the polaron band can be quanti- 
tatively described^ 3 , by the amplitude mode model 3 -^ of 
charge excitations in conjugated polymers. Features ii) 
or iii) are significantly weaker than the far-IR absorption 
i). We attribute the latter to the charge accumulation 
layer on the n-Si side of the gate dielectric. The validity 
of this assignment is supported by quantitative agree- 
ment of the far-IR data with results for heavily doped 
Si 3 -^ as will be described in detail below. 

We now turn to the analysis of the accumulation layer 
(AL) in n-Si. These results will help us to quantify the 
2-dimensional (2D) carrier density which is the same on 
both sides of the gate dielectric in the P3HT/insulator/n- 
Si structure. This task can be accomplished for exam- 
ple by fitting the low-energy absorption with the Drude- 
Lorentz (DL) model, in which the complex dielectric 
function s(lo) has the form: 



J pD 



iT 'aw 



iYuj 



(27) 



The first term in Eq. (|27|) stands for the Drude free car- 
rier response with the relaxation rate T d and the plasma 
frequency ujpp, given by Wp D = 47r ^ e , where Nf is the 
2D density of free carriers and mb = 0.26m e is the con- 
duction band mass of Si. The second term is used to 
account for small deviations of the dielectric function in 
mid-IR from the Drude form. The last term Soo is the 
high frequency dielectric constant. Two scenarios are the 
most likely origins of the second term: frequency depen- 
dence of the relaxation rate of free carriers and/or inter- 
band transitions from donor impurities to the conduction 



band of Si at loq. The interband transitions are charac- 
terized by the relaxation rate T and the plasma frequency 
lop given by cop = 47rjVj !f e , where Nm is the 2D charge 
density and m* ~ mt> is the renormalized mass. It is ad- 
vantageous to carry out the analysis using the "effective 
2D" carrier density in Eq. (|2~7| in order to avoid any as- 
sumptions about the thickness of the accumulation layer 
in n-Si. The total 2D carrier density A^d = N{+Nm rep- 
resents a lower limit of the density of the voltage- induced 
charges. 

In order to evaluate the consistency of the DL model 
with the data, we first obtained optical constants for 
each of the constituting layers of our devices as described 
in Ref. [34| and then considered the transmission of the 
multi-layered sample: P3HT/Si0 2 /2D-AL/n-Si. Results 
shown in the top panel of Fig. [3] suggest that the ex- 
perimental spectra can be adequately described by mod- 
eling the AL response using Eq. (|27[) with parameters 
presented in Table [TTJ Here we assumed that voltage 
induced changes of P3HT in far-IR are negligible small 
in accord with dc measurements^ The parameters we 
obtained for the response of AL are in good agreement 
with the optical spectra of metallic phosphorus doped 
Si^ The voltage-induced absorption spectra in the mid- 
IR frequencies for the P3HT layer can be obtained by 
subtracting the absorption due to the accumulation layer 
in n-Si from the absorption spectra of the device, which 
are displayed in the bottom panel of Fig. [3] 

In Fig. |4] we plot the 2D carrier density of charges A^d 
injected in P3HT polymers obtained from the previous 
analysis as a function of the gate voltage Vgs ■ The car- 
rier densities in the polymers reaches 10 13 cm" 2 at the 
highest voltage. Therefore, we have shown that the car- 
rier density in P3HT in our organic FETs at the highest 
gate biases is about A^ 3 d ~ 10 20 cm" 3 corresponding to 
approximately one charge in every 10 thiophene rings, 
i.e., 10%. 



V. DISCUSSION 
A. Comparison with the IR spectra 

Let us now compare the theoretical predictions with 
the IR spectroscopic results and discuss their implica- 
tions. In Sec. IIII M we have estimated the lowest optical 
absorption frequency to be 0.16eV for polarons and 0.46 
eV for bipolarons. In our IR spectroscopic studies, the lo- 
calized absorption band due to charge injection is found 
at about 0.45 eV. The frequency is appreciably higher 
than the predicted polaron absorption, while it matches 
reasonably well with the predicted bipolaron absorption. 
The result seems to suggest that the localized band is due 
to bipolaron excitations. On the other hand, polaronic 
energies reported in previous work seem to vary appre- 
ciably suggesting a dependence on the sample qualities 
and experimental platforms. Indeed, very broad range 
of polaron energies ranging from 0.06eV to 0.35eV has 
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FIG. 4: The 2D carrier density of charges and the correspond- 
ing doping level y (solid triangles, solid line is to guide the eye) 
induced in the P3HT polymer as a function of the gate volt- 
age determined as described in the text. The density reaches 
10 13 cm -2 (equivalent to 10%) at the highest Vgs applied. 
The horizontal lines mark the different ranges of the critical 
dopant concentration at which metallic states can be expected 
predicted from the theoretical models in Sec. IIII BfHTl Dl 

been reported in different P3HT systems4 ' 39 i 40 Consider- 
ing that the simple model followed in this paper assumes 
idealized ID infinite chains, there is a need to discuss in 
some details the disorder effect on the localized bands 
and the optical transition energies. This discussion is 
presented in the subsequent subsection. 

The calculations in Sec. IIIII show that a relatively low 
doping concentration of about 5% is expected to close 
the gap between the polaron bands and the continuum 
states, leading to a possible second-order transition in 
P3HT. In our actual experiments, with a gate voltage of 
— 100V applied, the charge injection of the P3HT poly- 
mers achieves a doping concentration as high as 10% (see 
Fig. 2]) , exceeding the threshold over which the polaron- 
related transition occurs. However, we do not have any 
direct evidence of metallic states, for example, through 
reduction of the oscillator strength of the IRAV modes 
near this doping range.— The insulating state at the dop- 
ing level of 5% is thus consistent with our bipolaron in- 
terpretation of the spectra. 

On the other hand, our results indicate that metallic 
states in P3HT can be expected at a doping concentra- 
tion of about 10% due to either a first-order transition 
to the metallic polaron lattice or a possible second-order 
transition by closing the gap between the bipolaron band 
and the valence band. Both thresholds are close to the 
highest doping achieved in our charge-injected polymers. 
In the following we suggest two issues that can lead to 
the absence of direct evidence of metallic signatures in 
the present experiment. 

First, recent experiments in P3HT FET suggested that 
at low temperature, a DC field is needed in order to over- 
come the tunneling barriers between the occupied local- 
ized states which can suppress conductivity of polymers 



and result in insulating properties^ Therefore while it 
is possible that the present P3HT system might have 
achieved a metallic carrier density, it is not clear that 
one can observe metallic conductivity, as found recently 
that only a small faction of the field-induced charge con- 
tributes to metallic conduction while most charges occu- 
pying the localized states^ Second, in the FET device of 
interest here, the absorption in the accumulation layer on 
the n-Si side of the gate dielectric dominates the observed 
spectrum in the far-IR frequencies, because of its much 
higher mobility compared to that of P3HT. As a result, 
it is difficult to directly subtract the Drude absorption 
due to the n-Si from the absorption of the device so that 
the Drude absorption and thus the metallic signatures of 
the polymer at the far-IR frequency can be extracted. 



B. Effect of disorder 

In quasi-lD conducting polymers, disorder is expected 
to have a particularly important role affecting both trans- 
port and other phenomena compared to systems with 
higher dimensions. Indeed, previous work has suggested 
that insulator metallic transitions in trans-PA can result 
from disorder in the polymers Here we discuss the sen- 
sitivity of our predictions to disorder of the polymers in 
terms of the confinement parameter 7 = y^-, which is 
directly related to the energies of the localized states. 

In the presence of a higher degree of disorder, we 
expect that the degree of the conformational nondc- 
generacy would be lowered due to the configurational 
randomness. Hence the external gap parameter A e 
is lower compared to that of more regio-regular poly- 
mers. More importantly, disorder affects the confine- 
ment parameters through the electron-phonon couplings, 
the strength strongly influences the electronic structure 
of organic polymers^ In organic semiconductors, the 
electron-phonon coupling constant is related to the con- 
formational disorder through the degree of the chain con- 
jugation. The coupling decreases with the conjugation 
length,— i.e., the length of the polymer segments sepa- 
rated by defects. Hence, the electron-phonon coupling 
increases with disorder in the conjugated polymers. On 
the basis of these arguments, one can anticipate that in 
the presence of disorder, the confinement parameters 7 
would be smaller than that in polymers with higher regio- 
regularity. 

In Sec. IIII Bl IHTDI we have carried out calculations of 
the critical densities in P3HT using 7 = 0.4. Given the 
trend that 7 decreases with increasing disorder, it is in- 
structive to recalculate the results at values of < 7 < 
0.4. To this end, we have recalculated the values using 
7 = 0.01 and obtained u = 0.75eV and ivo = O.leV 
that correspond to the lowest optical absorption ener- 
gies of 0.3eV and 0.95eV for the polarons and bipo- 
larons, respectively. Note that the two bipolarons bands 
merge into a single mid-gap state at 7 = 0. The results 
show that both the polaron and the bipolaron absorp- 
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tions blueshift with increasing disorder of the polymers. 
However, even at the smallest value of 7, it appears that 
the localized absorption in our experiment is still con- 
siderably higher compared to the predicted polaron ab- 
sorption, thus supporting our assignment of the localized 
band as due to bipolarons, rather than polarons. 

Extending this analysis to the larger end of 7, on the 
other hand, we find the lowest polaron absorption ener- 
gies reduced significantly to less than O.leV when 7 = 1. 
This value is consistent with the recent photoinduced ab- 
sorption measurement by Osterbacka et at of highly 
regio-regular P3HT polymers^ 3 - 9 - where the polaron ab- 
sorptions are found to superpose with the IRAV modes at 
very low frequencies less than O.leV. Similar absorption 
features were also revealed by Brown et al. using charge 
modulation spectroscopy in P3HT based MOS diode, al- 
beit assignment of the polaron absorptions were given 
somewhat differently^ It is important to note that a dis- 
tinction exists between the RR P3HT in our experiment 
and those in the earlier work ) 20 ' 39 despite both exhibit 
high mobilities. In our experiment, the charge accumu- 
lation layer is restricted to the polymer/gate dielectric 
interface within 1-2 nm of the 4-6 nm-thick P3HT films 
due to the bottom contact FET structure, the interfacial 
disorder— would play a much more important role than 
in the "bulk" form of P3HT in the earlier experiments. 
We thus suggest that the difference between the absorp- 
tion peak we observed and those from the earlier work 
might be related to the extent of the interfacial disorder 
present in the samples. 

We have also recalculated the doping thresholds at 
< 7 < 0.4. This leads to a range of the critical doping 
values as shown in Fig.HJ For the first two cases discussed 
in Sec. IIII Bl and IIII CI the thresholds increase slightly 
with decreasing 7, i.e. with increasing disorder. For the 
bipolaron gap closure discussed in Section UlI D[ the value 
of 7 is lower bounded by the value at which the two bipo- 
laron bands in the gap touch due to the band broadening, 
which corresponds to 7 = 0.3. In other words, below this 
value, the bipolaron gap closure condition can not be sat- 
isfied. In the range of 0.3 < 7 < 0.4 , the critical doping 
concentration is nearly unchanged. We have indicated 
the latter threshold by a single line in Fig. [4] Given 
the fact that these ranges appear relatively insensitive 
to the structural disorder described by 7, together with 
the assignment of the carries to bipolarons, we may con- 
clude that the maximum carried densities achieved in the 
present experiments are in close proximity to the regime 
where a first-order transition from bipolaron lattice to 
metallic polaron lattice occurs. 

We conclude this section by stressing that due to con- 
tact and interface effects, we expect the precise value of 
the critical density for insulator-metal transition to be 
larger than the one we have predicted above. In this re- 
spect, we compare our results to very recent experimental 
data obtained by Panzer et al. 15 and Dhoot et alM The 
first authors have reported field-induced charge injection 
as high as 10 22 cm~ 3 in P3HT FETs with polyelectrolyte 



gate dielectric. At such densities - higher than our pre- 
dicted minimal density - P3HT indeed shows metallic 
transport properties. The second group2i finds strong 
evidence of metallic behavior in a polymer which is struc- 
turally very similar to P3HT. This metallic behavior is 
again observed at carrier densities larger than the ones 
we predict. These results therefore corroborate our pre- 
diction that a critical carrier density exists beyond which 
metal-insulator transition can be expected in P3HT. 



C. Limitation of the model 

We briefly comment on the limitations of our theoret- 
ical analysis. The BK model used in this paper does not 
directly account for the electron-electron interactions. In 
spite of ongoing efforts to incorporate electron-electron 
interactions in the models of conjugated polymers^ it 
is currently far from clear how that has to be done. 
Nonetheless, it is worth noting that to some extent, 
electron-electron interactions are effectively included in 
our calculations through the use of experimental values 
for the parameters such as the band gap, the character- 
istic correlation length, and the confinement parameters. 
We expect that explicit inclusion of the interactions in the 
present treatment may modify somewhat the critical con- 
centrations predicted in the analysis. On the other hand, 
the theoretically calculated polaron/bipolaron energies 
and the experimental values obtained from our P3HT 
FET are in reasonable agreement. This agreement and 
the insensitivity of the results to the doping concentra- 
tion, and hence to the strength of the electron-electron 
interactions suggest that the electron-lattice interactions 
are still the dominating effect in P3HT polymers. The 
electron-phonon coupling model we use thus captures the 
essential physics of localized excitations. 

The BK model also does not account for the electric 
field. While the effect of the electric field is important for 
understanding polaron dynamics, we expect that it does 
not play an important role in regioregular P3HT-based 
devices considered in this paper. In the P3HT structures 
we have studied, the P3HT chains form lamellae normal 
to the substrate with the chains oriented parallel to the 
substrate surface. 8 As the applied gate field is perpendic- 
ular to the film surface, we expect that the P3HT chains 
are not subject to the gate field in the parallel direction 
and thus the electrostatic potential along each chain is 
constant. 



VI. CONCLUSION 

In this paper, we have analyzed the conditions to ob- 
tain metallic behavior in charge-injected P3HT polymers 
through a first-order phase transition between bipolaron 
and polaron lattices or gap closure between txshe bipo- 
lar/polaron bands and the continuum bands. IR spec- 
troscopic measurements in P3HT thin film FET show 
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that the highest charge injection level achieved in the 
present experiments is in proximity to the predicted dop- 
ing thresholds for the first case, indicating that metallic 
behavior can be anticipated. Optical absorption frequen- 
cies are predicted for both polarons and bipolarons. The 
frequencies associated to bipolarons are in good agree- 
ment with the absorption peak observed. Both conclu- 
sions remain true even when disorder is taken into ac- 
count. We therefore conclude that the localized excita- 
tions we observe are likely to be due to bipolarons and 



a first-order transition from the bipolaron lattice to the 
polaron lattice can be expected. 
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